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ABSTRACT 

We derive crystal braneworld solutions, comprising of intersecting families of parallel n+2- 
branes in a 4 + n-dimensional AdS space. Each family consists of alternating positive and 
negative tension branes. In the simplest case of exactly orthogonal families, there arise 
Oh. different crystals with unbroken AD Poincare invariance on the intersections, where our 
r-| ! world can reside. A crystal can be finite along some direction, either because that direction 
is compact, or because it ends on a segment of AdS bulk, or infinite, where the branes 
continue forever. If the crystal is interlaced by connected 3-branes directed both along 
the intersections and orthogonal to them, it can be viewed as an example of a Manyfold 
universe proposed recently by Arkani-Hamed, Dimopoulos, Dvali and the author. There 
are new ways for generating hierarchies, since the bulk volume of the crystal and the 
lattice spacing affect the 4D Planck mass. The low energy physics is sensitive to the 
boundary conditions in the bulk, and has to satisfy the same constraints discussed in 
the Manyfold universe. Phenomenological considerations favor either finite crystals, or 
crystals which are infinite but have broken translational invariance in the bulk. The most 
distinctive signature of the bulk structure is that the bulk gravitons are Bloch waves, with 
a band spectrum, which we explicitly construct in the case of a 5-dimensional theory. 
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Recently a remarkable observation by Arkani-Hamed, Dimopoulos and Dvali that there 
may exist additional sub-millimeter spatial dimensions [|1| has generated tremendous in- 
terest. In such theories, the usual Standard Model degrees of freedom are localized to a 
3-brane which is embedded in a higher- dimensional bulk, and thus at low energies are in- 
different to the extra dimensions. On the other hand, gravity and, typically, other weakly 
coupled fields live in the bulk, probing equally all spatial dimensions. The observed weak- 
ness of gravity in four dimensions is generated naturally, because by spreading through 
all spatial dimensions gravity becomes softer, as can be seen immediately from Gauss law. 
Hence the hierarchy between the Planck scale, Mpi ~ lO^^GeV^ and the electroweak scale, 
'^EVF ~ TeV and its radiative stability, can be naturally explained if the size of the n 
internal dimensions r is large: Mp^ = M^+^r" ^ ~ m\-^ 

Studies of phenomenological constraints in |^ have confirmed that such models are 
consistent with observations, with the unification scale as low as few x 10 TeV . There 
has subsequently been much interest in models with large extra dimensions String 
theory may give rise to such models, with the radius of compactification only a few orders 
of magnitude above the Planck scale, such as in Hofava-Witten p!0| , |Tl|, intermediate sizes 
~ (TeV)-^ PI g, and even ~ {Mpi/ Mf''' M-^ < mm | 



The concept of localization of matter to the branes plays a key role in masking the 
large extra dimensions from observation at low energies. In fact, the universe as a domain 



wall in a non-compact space has been considered before |T^, |T6|, but there gravity was 
higher-dimensional at all scales. But if gravity itself is localized to the 3-brane, then 
AD gravity can be reproduced at low energies even if the extra dimensions are much 



larger than a mm. A very interesting recent construction by Randall and Sundrum [|T7 
provides an elegant example of this idea for the case of one extra dimension, which is 
linearly infinite, but has finite proper volume. Similar ideas have been explored also in 
| |18[ |. Subsequently it has been shown that gravity can be localized to intersections of 
n + 2 branes in 4 + n-dimensional spaces . More examples were found later , and 



various issues were considered. Since the Randall-Sundrum proposal [0 employs 

AdS bulk geometry, it is possible to view the model in the framework of the AdS/CFT 



correspondence p6|, Interpreting the fifth dimension as the scale of the AD theory 



p8[1 and the evolution along it as the holographic renormalization group [^, this model 
can be rephrased as a AD CFT theory coupled with gravity [Q, without invoking extra 
dimensions. The AdS/CFT correspondence has also been pursued in pT[ . 



On a different track, it is possible that our 3-brane is not alone in the bulk. In fact, 
the bulk may be populated by many other branes, if we wish to understand the weak 
breaking of symmetries by shining in the bulk P] . Hence, the universe may contain many 
nearby branes, which are dark because the light travels only along them, but are felt 
by their gravity ||32|. This Manyfold universe has been considered recently by Arkani- 



Hamed, Dimopoulos, Dvali and the author [^, and it has been shown that it can be 
consistent with the observations, while giving rise to many new phenomena: dark matter 
candidates which are dissipationful, dark stars, hybrids, neutrino mixing, SUSY breaking 
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et cetera. Some solutions representing multi-brane configurations in 5D spacetimes have 
been constructed and some aspects of brane networks considered [0, p5[] . 



In tfiis paper, we will construct explicit solutions of Einstein's equations which repre- 
sent intersecting families of parallel n + 2-branes in 4 + n-dimensional AdS spacetimes. 
Each family will contain branes of alternating tension. While we are a little reluctant 
to treat the negative tension branes on equal footing with the positive tension ones, we 
will ignore this when considering the solutions because of considerable technical simpli- 
fication. If the bulk cosmological constant is replaced by a potential, there may exist 
solutions which involve only positive tension branes, while sharing some generic prop- 
erties with our solutions. Alternatively, a first-principles construction of branes with 
negative tension may yet ensue. We will show that there arise several different crystal 
solutions, all of which support the usual AD gravity. However phenomenological consider- 
ations strongly favor finite crystals or crystals which are infinite but have strongly broken 
translational symmetry, where our world arises as a lattice defect. If a crystal is inter- 
laced by a grid of 3-branes, some along the intersections and others stretching between 
them, the low energy physics of such a solution is identical to that of a Manyfold universe 



32]. Such a crystal would have to be subject to the same astrophysical constraints as 



any generic Manyfold . We will further show that the spectrum of bulk gravitons has 



band structure, with a gap at zero energy, because the bulk gravitons are Bloch waves, 
like electrons in ordinary crystals. This may allow the string scale in a 5D theory to be 
as low as M^: ~ lOOTeV^, at the cost of having Af ~ 10^^ branes in the 5D AdS space- 
time. However, the reward consists of reproducing the AD Planck scale Mpi ~ lO^^GeV, 
the electroweak scale tjiew ~ ITeV and suppressing the corrections to Newton's law at 
distances > 1mm. It may be possible to reduce the number of required branes by going 
to more than one extra dimensions. 

We begin by deriving the solution describing the crystal braneworlds. Consider an 
array of n orthogonal families of n + 2-spatial dimensional branes in a 4 + n-dimensional 
spacetime, with a bulk cosmological constant A. For simplicity we take the branes in each 
family to have identical tension ak- The field equations can be derived from the action 

S= f rf'+"a;y^(-^ + a) - E E / ^'"""a^y^afc. (1) 

Here n^^n — 8^/^"'*'^ 5 where is the fundamental scale of the theory. In the brane 
actions, refers to the induced metric on the brane. Note that the measure of integration 
differs between each brane, and between the branes and the bulk. This will be reflected 
in the field equations, where ratios weigh the (5-function sources. The gravitational 

field equations are 



G% = kI^^A5\ - E ^^UK4+n^fc E(-l)'^(^' - Jh)die.g{l, 1, 1, 1, 1, 0,, 1) (2) 

where the coordinates z'' parameterize the extra dimensions. In general, these equations 
should be supplemented by the equations of motion for the branes, which however are 
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solved by placing branes at fixed locations parallel to the AdS boundary. The issue 
of stability would then have to be addressed, and is equivalent to providing the brane 
G oldston e bosons with mass terms. We will not address this in detail here. The ratios 

'^^^"''^^ which appear in (H) reduce to ^Jg^ for simple diagonal metrics, however in 
general they cannot be gauged away. 

Deriving the solutions is straightforward ||TP[. Away from the branes, the bulk geom- 
etry consists of patches of the 4 + n-dimensional AdS space. On the branes, there are 
additional (5-function singularities, which cause discontinuities in the derivatives of the 
metric. To find the global solution, valid on the complete manifold, we should seek the 
solution in the form 

dsl^, = [ii^.dx^dx" + (^'^')') (3) 

k=l 

and adjust the boundary conditions for Qr^ such that it solves Einstein's equations (0), 
and between the branes reproduces the metric of AdS bulk. It is convenient to recast 
the Einstein's equations in the conformal frame Qah = ^~^gah- We use the standard 
relation for 4 + n spacetime dimensions 



Gab = Gab + {n + 2)(yalogQVblogQ-VaVblogQ) 

+ {n + 2)~ga, (V'logn + "^[Vlognf] . (4) 



Since the conformal metric is flat. Gab = 0. Then substituting (^) into (0) yields two 
equations for VL~^: 



72n-l _ 



{vn~')' = 2 , ^'+;^ (5) 

^ ' in + 2){n + 3) ^ ' 



Note that in deriving the first equation it was important to use = ^, which 

ensures that the brane tensions are constant. 

It is now clear from the first of equations that Vl^^ is a linear combination of Green's 
functions for each family of parallel branes. Since the tensions 0"^ are the strengths of 
sources, they control the discontinuity of the derivatives of VL~^ on each brane. It is easy 

to see that the discontinuities must satisfy = -^i^^- Thus the discontinuity in 

the gradient of on each intersection of branes must satisfy (VVt^^Y = 4(^^|yT- But 
because the second of eqs. is defined globally, these two expressions must coincide, 
leading to 

4 2 _ 8(n + 2) 2 . ..^ 
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Hence the condition for the existence of the solutions is to fine tune the tensions of the 
branes to the bulk cosmological constant according to (^. Then as we see from the 
induced metric on the intersections is Minkowski, meaning that the conditions (^) amount 
to fine tuning the effective 4D cosmological constant on each intersection to zero. 

We can now solve the first of (||) as follows. Suppose that in each family there is A4 
parallel branes, which for reasons of simplicity we have taken to be equidistant, separated 
by Ik, as is seen in eqs. and (|^). The solution can be easily generalized when the 
distances between the branes vary within each family, as will become clear from the 
construction below. The eqs. (H) yield the solution for the conformal factor Q 

n-^ = Kj2^iz'') + i (7) 

k=l 

where K = {y/nL)~^, and the functions S{z^) satisfy 

J2 Q(^k 



k\ 



1. (8) 



These two equations follow from (||) by using linear superposition, eq. and the def- 
inition of the AdS radius. The integration constant in must be nonzero in order to 
excise the AdS boundary from the manifold, and is set to unity by gauge choice. As a 
result, the warp factor never diverges; otherwise there could not exist a normalizable AD 
graviton mode. 

It is easy to see that the function iS(z^) which solves eqs. must be the sawtooth 
function. We can write it as follows: 

_ I "^P^^ ~ " <z^ < 2plk; 

'^^^ ^ - - 2ph, 2ph <z'< {2p + 1%; 

The solutions differ globally depending on the boundary conditions for each In 
general there are four distinct types of boundary conditions: 

(i) There is a discrete infinity of branes along the k^^ axis, and the crystal is infinite 
in this direction; 

(a) There is A/fe = 2Nk + 1 branes, and N^. + 1 have positive tension and negative; 
the crystal ends on branes with positive tension in each direction along the z^ axis; outside 
of the crystal there is the near-horizon geometry of AdS space, with the horizon located 
at infinite proper spatial distance from each end; 

{iii) The crystal is semi-infinite along the z^ axis, with infinitely many branes, and 
boundary conditions satisfying (i) on one end and (ii) on the other; the brane at one end 
must have positive tension; 
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(iv) The crystal is finite along z'^, with an even number A4 = 2Nk of distinct branes 
along the half of them with positive and a half with negative tension; the 

solution can be viewed as case (ii) but with the 0**^ and 2Nk + 1*^ branes identified; this 
means, that the axis is compactified on a circle. 

By linear superposition, the global solution for the crystal is a combination of any of 
the four possibilities realized along individual axes. For the case (i), the solution for the 
sawtooth functions is adequately represented by an infinite array of the form (^). In the 
other three cases, by choosing the origin of the coordinate system in the bulk to lie at the 
end brane with positive tension, we can rewrite the solution in a more compact form: 

S{z'^) = 2Y:i-iy9{z'~jlk){z' -3k) - z\ (10) 
i=o 

where Mk = 2Nk + 1 in the case (ii), Afk ^ cxd in the case (in), and Mk = 2Nk in the case 
{iv). Therefore, in terms of these functions, the metric can be written down as 

1 " 



This is the solution for the crystal Manyfold universe. By the conditions (|^), all the 
vertices of the crystal lattice correspond to the 4D Minkowski spacetimes, and can ac- 
commodate 3-branes which localize the Standard Model particles. Therefore, our world 
could reside at any one of the vertices, leading to different low energy theory, as long 
as there exists a massless AD graviton mode. The 3-branes lying on the other vertices 
would behave as mirror worlds, and can be connected to our own by additional 3-branes 
embedded in the intersecting n + 2-branes. Hence a crystal interlaced by 3-branes would 
give rise to a Manyfold universe structure of p2|. 



It is straightforward to verify that the theory defined by the action (|T]) in general admits 
AD massless graviton modes on the backgrounds ([TT|) . Hence at large distances these 
modes will give leading contributions to the gravitational interactions between particles 
localized to the intersections, reproducing the usual Newton's law and masking the extra 
dimensions. The wave function of these modes is 

Vo(a;^i')~^^('-")/'(^)e^'^''^^ (12) 

as we will see below in more detail. Using (|l]) and (|12]), by ignoring the massive KK 
modes and reducing the theory (|l]) to four dimensions, we find that the 4D Planck mass 
is 

Mli = M:^+^f (13) 

J crystal space 

The integral has to be performed over the whole crystal space, which is composed of 
a number of patches of a different kind. Depending on the type, there may be semi- 
infinite patches comprised of the near-horizon segments of AdS bounded by branes of 
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positive tensions. Each such patch is identical to the patches considered in the case of a 
single intersection of branes, contributing (^n+iy.K^ ^*^"^ ^pi tUll- However, since now 
there are also patches of geometry which are completely bounded by branes, there arise 
additional contributions to Mpi. Different patches, or cells, can be classified by counting 
how many branes they are bounded by, such that the integrations in those directions 
would be bounded between and 1^. It is clear from (|^) and (y) that all topologically 
identical segments give the same contribution to the mass integral. Each distinct integral 
needs to be weighed by the symmetry factor, which counts the total number of cells of a 
given class in the crystal. The integrals I = /^gn d"'z fi""*"^ in general assume the form 

where p is the number of directions which end on the AdS horizon. The integral can be 
evaluated straightforwardly, giving 



71— p 
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where Bj are combinations of {/p+i, of length j. To compute the symmetry factors, 
we break up the crystal into topologically inequivalent cells. These are distinguished by 
how many of their sides end on AdS horizon. Clearly, the cells inside the crystal only 
end on branes, and correspond to taking p = in ([T^). We can count the cells in each 
topological class as follows. In the case of ID crystals, we can divide the cells along the 
crystal into two types, A and B. The type A are the cells inside the crystal, which end on 
branes on both sides, and there is A4 of them. The type B cells end on the AdS horizon, 
and there is at most two of them. Then the crystal ( |lTl) can be represented formally as a 
direct product 

n 

Brane Crystal = \{ [NjA + QB] , (16) 
(g) i=i 

where the number of end branes is C = 0, 1 or 2 depending on the boundary conditions. 
It is then obvious that the number of elementary cells which end on the AdS horizon in 
p directions is given by the coefficient of B^ in the expansion of (pIBI) : 

cEn[^], (17) 



p 



c 



where Cp is any set of n — p distinct elements from the set {l,...,n}, and [r] denotes 
rounding off the number r to the nearest smaller integer. Specifically, the number of 
cells inside the crystal, with p = is z/q = C" 11^=1 ["^l- Note that the total number of 
branes along k^^ direction is set by the boundary conditions {i)-{iv). Hence we obtain the 
following formula for the 4D Planck mass: 

/■nnfn+2 n sT n-p . 

*'"=(^Tw^^nrai:(-i''T7r-7f^^^. (18) 
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where Bj are combinations of length j of elements of {lc{i), ■■■Jc{n-p)}, and c{k) G Cp. In 
the special case of finite crystals with cubic global symmetry and two end branes, there is 
the same number of branes in each direction, Afi = ... = Mn = A/", and the lattice spacing 
is the same in all directions, /i = ... = /„ = /. Therefore all Cp's of the same length 
coincide, the symmetry factors Vp are Vp = 2"[:|^]'^~p ^,^^j_^^, , and hence 

It is now obvious that for either class of finite crystals, (ii) or (iv), the formula (0) gives a 
finite value for the AD Planck mass, meaning that there is a 4D graviton. In effect, those 
crystal configurations effectively "compactify" the extra dimensions. In the case of infinite 
crystals (i) and {iii), although each individual contribution to ([T8|) is finite, formally the 
sum diverges since there is an infinite number of cells to sum over. However, the sum 
can be renormalized at the expense of loosing a one-to-one correspondence between the 
4:D and {4 + n)D Planck scales. Indeed, for infinite crystals Mpi ~ M^M'^^'^L^ oo, 
since A/" oo. Since the divergence comes from the bulk volume being infinite, the 
renormalization must be a bulk phenomenon. To renormalize Mpi, we recall that since 
by the fine-tuning condition (^) the bulk cosmological constant is related to the brane 
tensions, and since the brane tensions must be unaffected by a bulk renormalization, the 
AdS radius must also be unaffected. Hence, the only way to remove infinity from Mpi is 
to take M^, to be infinite, and renormalize it by Mbare A/'^"/'-"^^''M*. But this means 
that the resulting mass M^, need not be related to the scale of the brane physics, and 
hence this would seem to require another fine tuning to explain the hierarchy between 
Mpi and itiew- 

On the other hand, eq. ([T9|) shows that for finite crystals there are additional possi- 
bilities for generating hierarchies. This is because Mpi depends on four different scales: 
the fundamental scale M^, the AdS radius L = {\/nK)~^., the lattice spacing / and the 
total volume of the crystal ~ {Ml)"'. Hence it is easier to generate exponential hierarchies 
in large crystals, even if all the other scales are of the same order of magnitude. 

Before we move on to study the bulk gravitons, we should note that infinite and semi- 
infinite cases {€) and {Hi) are also disfavored phenomenologically, even if we perform a 
renormalization as outlined above. In an infinite crystal there is an infinite number of 
intersections. Discrete translational invariance of the crystal would require that the energy 
density of particles on each intersection is of the same order (which would remain true even 
if the symmetry is weakly broken). Because at long distances the gravitational force is 
dominated by the zero modes, an observer on any intersection would feel the gravitational 
field of particles localized to all the intersections in the crystal regardless of how far in 
the bulk they are. Thus the localized energy density everywhere would add up to give 
the total energy density controlling the cosmological evolution of all intersections. Hence 
to sum up to a finite value, the energy density on each intersection must be infinitesimal. 
This implies that all but an infinitesimal amount of the energy density of the Universe 
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we experience would have to be dark. Note that this remains true regardless of any 
renormalization of the scales, which may or may not be necessary (in a fashion similar to 
that for Mpi). The matter on distant intersections would still be dark, there would be 
infinitely many of them, and renormalization would not change the ratio of dark matter to 
visible matter. Since the amount of visible matter in our Universe is not infinitesimal, we 
see that at least one intersection in the crystal must have finite p, and most of the others 
cannot. Hence the discrete translational invariance of the crystal must be strongly broken. 
Our intersection is a defect of the lattice. It cannot deform the lattice too strongly, in order 
not to loose the AD graviton. But it does suggest that the graviton zero mode accumulates 
around the defect, which deforms the background more than most of the other vertices 
in the lattice. This can generate a peaked profile for the wave function of the zero mode, 
making other intersections weakly coupled. Such possibilities for generating asymmetric 
mirror worlds have been discussed in [^]. In practice, therefore, an infinite crystal would 
have to be physically hardly distinguishable from a finite crystal, with the exception that 
the strong coupling regime, discussed in ||T^, could be absent. A finite crystal would 
have to be consistent with observational constraints much like the Manyfold universes 
in theories with sub-millimeter internal dimensions Thus, an infinite crystal with 



unbroken discrete translational invariance may be only a mathematical idealization. We 
will retain it, however, since it is quite useful for the study of bulk gravitons, to which we 
now turn. 

The field equation for linear perturbations around the background solution can be 
obtained by expanding the 4Z) part of the metric according to g^^ = g^y + h^y, where h^i, 
is in the transverse traceless gauge with respect to the background V ^h^^y = h^^^ = 0, 
and expanding the Einstein's equations to linear order in hfj_y. The perturbations satisfy 

6R^y = ^R\h^x + ^R\hxy. (20) 

After evaluating SRfj^y, substituting the background crystal metric (PD, and defining the 
wave function by 

V = (21) 
where in general denotes a complex function, which arises as a linear combination of + 
and X polarizations, we find 

1^ T 1^2 T ra(n + 2)(n + 4) ^ 



2 2 8{EkS{z^) + V^L) 



2 



^ 2(E.;(")lv^L) SB-im^--A)^^0. (22) 

In this equation, the 4D D 'Alembertian is solved by the expansion into plane waves, 
□4\E' = m^\E', where m is the KK mass. The resulting equation is 

v^* + H - + ,^ rn t U-m^' - A))* = 0. 

(23) 
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This equation is a static Schrodinger equation for an electron in a periodic potential. It 
is now clear that the lowest energy eigenfunction of (0), which corresponds to the 4D 
graviton zero mode, is ~ const with m? = 0, which by ( pi]) reproduces precisely ([T2|) , 
as we claimed. For > 0, the analogy with the electron in a crystal shows that the 



solutions of (^31) are Bloch waves, satisfying ^q{z) = exp{iq- z)v^{z), where by Floquet's 
theorem are periodic functions under the action of the discrete translation group. 
The energy spectrum has a band structure due to resonant scattering. The bands are 
continuous in infinite crystals, and discrete in finite ones. In addition, there is a mass 
gap, separating the lowest-lying excited state from the ground state by a finite energy. 
In the language of our crystal Manyfold, this means that the lightest KK graviton must 
have a finite mass, and that the KK modes come in bands. As a result, their influence on 
the low energy physics on the intersection will be significantly suppressed when compared 
to a single brane of ||T^ or a single intersection of [|T9|. This is true not only for the 
production of bulk gravitons by the collisions of the Standard Model particles, but also 
for the corrections to Newton's law of gravity. Due to the gap, there will be an exponential 
suppression of the higher-order corrections to the inverse square law. Mass gap was also 
and in 



discussed in 23 



35 



To be more concrete, we focus on an infinite crystal with only one internal dimension. 
Finite crystals will behave similarly, when they consist of many branes, except that the 
allowed bands will be discrete rather than continuous. We expect that the qualitative 
properties of crystals with more than one internal dimension will remain similar. In the 
case of one extra dimension, the Schrodinger equation (|23|) reduces to 



15 



4(5(z) +L)2 



S(z) + L 



(24) 



By the analogy with electrons in a crystal, e.g. with the Kronig- Penney model, we need to 
solve this equation inside two adjacent elementary cells in the crystal. By the periodicity 
of (0), we can choose the elementary cells < z < 21 and 21 < z < 41. The 5-function 
potentials in ( p4D can be recast as the boundary conditions on the derivatives of \Ef at the 
vertices: ^'(/+) - ^'(/_) = and ^'(2/+) - ^'(2/_) = -f ^(2/). Moreover, the 

wave function at each vertex must satisfy {vertex-^) = {vertex-) in order to conserve 
probability. The differential equation in the first elementary cell is 



4{z+L)2 ^' 
15 

4(2/-2+L)2 ' 



if < z < /; 
I <z <2l. 



(25) 



These are Bessel equations. Using the boundary conditions above, the solution is 



m{z + L) AH^m{z + L)) + BH^imiz + L)) 



m{2l - z + L) 



^B]Htm{2l 



A + 



K^2 



^B)H2{m{2l 



-z + L)) 
z + L)) 



< ^ < /; 

1 <z<2l. 

(26) 
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Here are the Hankel functions, which generahze the plane waves for the case of the 
Bessel equation. They should be used instead of the real Bessel functions, since they 
are needed to encode the momentum by phase rotation. It is rather amusing that both 
classical 4D graviton polarizations are needed to mimic the quantum-mechanical phase. 
Bloch waves physically correspond to rotation of the polarization plane of bulk gravitons 
in the crystal. In evaluating the transmission and reflection coefficients in (EBI) on the 



5-function, we have used a recursive relation = Ht — —Hf. Then we can evaluate 

' aw w ^ 

the functions and derivatives at l±, and using the boundary conditions eliminate C and 
D in favor of A and -B, and have the function in each elementary cell be determined by 
two (complex) constants, just like in the Kronig- Penney model. 

The solution in the adjacent cell can be obtained similarly. It is given by 

'AHtm{z -21 + L)) + BH^im^z - 2/ + L))] , 2/ < z < 3/; 



m{z -21 + L) 



m{Al - z + L) 



3/ < 2 < 4/. 



+ [tt§^A + 'i^^i^^B)H^{m{Al -z + L)) 



(27) 

Now we need to relate A, B, A, B. The boundary conditions which they must satisfy are 

*(2/) =exp(2zg/)^(0), 
^(4/) =exp(2zg/)^(2/), 
^{21+) = vl/(2/_), 

*'(2/+) - ^'(2/_) = -|^(2/), (28) 

where the former two come from the Floquet's theorem for the Bloch waves and the 
periodicity of the braneworld crystal, and the latter two come from the continuity of 
and the jump of at z = 21. The parameter q represents the wave vector of the Bloch 
wave, which is to be determined. The first three conditions yield 

^ _ {hih^ + h];h^)h^ - {h{h^ - h^h2) exp{2iql)h^ - 2h'[hth2 
{hlh^ + hih2)h2 + {hih2 — hih2) exp{2iql)h2 — 2/ij~/;,2 ^2" 

^ _ {hih2 + h^h^)h^ - jh^h^ - /i^^/ia ) exp{2iql)hi - 2h'[h^h2 
{hih2 + hih2)h2 + {hih2 — /i]^/i2 ) exp{2iql)h2 — 2hih2h2 

A = exp{2iql)A, (29) 

where we have defined the numbers = H^{m{l + L)) and = H^{mL). The last 
condition will give an equation relating q and m, because the normalization constant 
A factors out: it is determined by the overall normalization of the wave function. The 
equation for q is the band equation. A straightforward computation gives 

^ (j2ni + jln2)(j2^^l + jiri2) - jij2(jij2 + 3nin2) - nin2(3jij2 + nin2) 

cos(<gj = . (3(Jj 

2(j2ni - iin2){i2ni -jih2) 
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In this equation, we have decomposed the Hankel functions into the Bessel functions 
Jfc and Nk, = Re H'^^Nj. = Im H^, with the definitions jk = JA:(m(/ + L)), Uk = 
Nk{m{l + L)), jk = Jk{fnL) and = Nk{mL). To see that there is a gap in the spectrum 
aheady at zero momentum, we can substitute g = in (|30|). The resuhing condition for 
the lowest-lying mass m is the gap equation: 

jlj2jlj2 + nin2nih2 + J2j2^l5^1 + jljl^2n2 = 3jij2^1^2 + 3jlj2^l5^2 - 3jlj2^1?^2 - 3jij2nin2. 

(31) 

It is evident that this equation provides a non-trivial constraint on m. Finding the lowest 
positive value of m which solves this equation is a straightforward, albeit tedious task. 
However, for I > L we can readily obtain the correct order of magnitude of the gap. Since 
the Bessel functions in ( pT]) are evaluated at m{l + L) and mL, the only two candidates 
for mass scales for mgap are and (/ -|- L)~^. On the other hand, in the limit / — >■ cx), 
the crystal solution reduces to the single brane solution of |T^, where the bulk gravitons 
do not have a gap. This excludes L^^ . Hence the mass gap must be 

^gap = (32) 

With this at hand, it is easy to estimate the corrections to Newton's law in AD induced 
by bulk gravitons. By the form of the wave functions (p6D, the AD Newtonian potential 
for particles localized to a 3-brane at the intersection of 4-branes with positive tension 
("PlanckQ brane"), with the corrections from the bulk gravitons, is 

V„ = -G„^^i^(l+ar ^^e-). (33) 

r ^ Jmgap K K ' 

which upon integration gives 

V, ^ -G„^(l + a-4!-e-'"/<'«>), (34) 

where a and b are constants of order unity. Note that because of the gap, the corrections 
are 0(r~^), in contrast to (9(r~^) in []T^. On the other hand, using eq. (P!T|) when / > L, 
the potential for particles localized to a 3-brane at the intersection of branes with negative 
tension ( "Tel^ brane" ) is, roughly, 

yN = -G.^{l + al^^^ — {-)e ), (35) 

which gives 

l-.^-G.^fl + a-^e--/'), (36) 

and again a and b are (different) constants of order unity. Clearly the corrections to New- 
ton's law for particles on the Planck brane can remain small at sub-millimeter distances 



^The nickname "Planck" refers to the 5Z? Planck scale Af*, not the 4Z? Mpi. 
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even if L ^ /p^ ~ 10~'^^mm, despite the fact that there is only five spacetime dimensions. 
However, the corrections become more significant on the TeV brane, as is clear from (|36|). 
But if we are to solve the hierarchy problem while suppressing the corrections to AD 
gravity to satisfy the observable bounds, our world should be a TeV brane. Therefore, 
the existing observational constraints yield < 1mm ■ L^, i.e. / < 1mm ■ ( j)^ < 1mm. 
Further, if / > L, eq. (^) gives Mp^ ~ M^LAf, and so the AD Planck mass depends 
on the size of the crystal and the lattice spacing. The fundamental scale should satisfy 
{■ TeV. Hence M|,; ~ {TeVf A/", and using / ~ eV ^ <C 1mm, for which 
L ~ 10^^ ensuring the validity of the supergravity description of the solution ([TI|), 
gives Af ~ 10^^. With this choice of numbers, the Standard Model on the TeV branes 
has scale mEw ~ ITeV, and its mirror on the Planck branes has scale M^, ~ lOOTeV. 
While the number of branes required is clearly very large, being equal to i.e. to the 
hierarchy between the scales, it is smaller by sixteen orders of magnitude than the max- 
imal number of branes which fit inside the flat extra dimensions Q . Note also that the 
size of the crystal is A/"/ ~ lO^^mm, with the choice ~ lOOTeV, which is four orders 
of magnitude smaller than the size of a fiat compact dimension required to generate the 
same ratio . The required number of branes could be reduced in the case of more 

extra dimensions, but the details are beyond the scope of the present article. 

It is straightforward to generalize this discussion when the lattice spacing varies within 
each family of parallel branes. In the extreme case, the configuration may be more 
similar to few separate crystals residing at different places in the AdS bulk. The effective 
scale would change from crystal to crystal, and in the infrared limit the low energy AD 
physics on distant crystals would be similar to that on a 3-brane far from an intersection 
Toll , or far from the Planck brane in 5D |22|, or to a Manyfold with bulk profiles 



With several crystals in the bulk, it is easier to generate hierarchy with fewer branes in 
each crystal, as is straightforward to verify. In fact, separated "small crystals" would 
correspond precisely to the cases of distant branes discussed in |jl9|, |22|, where hierarchy 
is generated by gravitational shining. 

We close with several concluding remarks. We have constructed general solutions 
describing intersecting families of parallel n + 2-branes in 4 + n-dimensional AdS space. 
The families intersect at right angles, and consist of branes with both positive and negative 
tension. As a result, there arise several different types of crystals, both finite and infinite. 
All of them give rise to the usual AD gravity, but phenomenological considerations strongly 
favor either finite crystals or infinite ones with broken translational invariance, profiles in 
the bulk and our world as a lattice defect. If such solutions are interlaced by 3-branes, 



at low energies they behave the same as Manyfold universes of The periodicity of 
the crystal lattice implies that the bulk gravitons are Bloch waves, with a band spectrum 
and a gap at zero energy. We construct the spectrum explicitly in the case of one extra 
dimension, and find that in the extreme case it may be possible to have the string scale at 
~ lOOTeV^, if there is ~ 10^^ branes in the configuration. The crystal worlds raise many 
interesting questions: the issues of stability, early cosmology, etc. While the bulk gravitons 
behave as electrons in the usual crystal, one may wonder which degrees of freedom play 
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role of the conventional phonons. It seems plausible that they should be the Goldstone 
modes associated with the branes in the crystal, which may get excited by sufficiently 
strong external perturbations. A more precise consideration would be needed to test this. 
Further, it would be interesting to find the properties of higher-dimensional crystals, since 
they may require fewer branes to resolve the hierarchy problem while ensuring Newton's 
law of gravity in 4D at distances larger than a millimeter. We hope to return to some of 
these issues in the future. 
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